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Abstract. This article is intended to give an overview on a collection of results 
which have been obtained jointly with I.I. Ikromov, and in parts also with M. Kempc, 
and at the same time to give a kind of guided tour through the rather comprehensive 
proofs of the major results that I shall address. All of our work is highly influenced 
by the pioneering ideas developed by E.M. Stein. 

Consider a smooth hypersurface S in K 3 of finite type, and let dp = pda be a surface 
carried measure with smooth, compactly supported density p > with respect to the 
surface measure da. The problems that I shall address are the following ones: 

A. Find, if possible, optimal uniform decay estimates for the Fourier transform of 
the surface carried measure dfi. 

B. If we denote by A t the averaging operator A t f(x) := J s f(x — ty)dp(y), deter- 
mine for which exponents p the associated maximal operator 

Mf(x) := sup \ A t f(x)\ 
t>o 

is bounded on LP(R 3 ). 

C. Determine the range of exponents p for which a Fourier restriction estimate 

\f(x)\ 2 dp(x)Y /2 <C\\f\\ LP(R3) 
s ' 

holds true. The first problem is classical, and the other ones originate from seminal 
contributions by E.M. Stein. These problems are somewhat interrelated, and in 
particular it is well-known that the answers to B and C are strongly influenced by 
the answer to A. Moreover, due to the insights by V.I. Arnol'd and his school, in 
particular through work by A.N. Varchcnko and V.N. Karpushkin, we know that the 
answer to problem A can be given in terms of Newton polyhedra associated to the 
given hypersurface S, at least for analytic hypersurfaccs. That Newton polyhedra 
would play an important role also in other problems of harmonic analysis has been 
stressed later in pioneering work by D.H. Phong and E.M. Stein on oscillatory integral 
operators. 

In this article, I shall outline how (almost) complete answers to the above questions 
can be given again in terms of associated Newton polyhedra (for question B, at least 
if p > 2). 
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1. Introduction 

Let S be a smooth, finite type hypersurface in M. 3 with Riemannian surface measure 
da, and consider the compactly supported measure dp := pda on S, where < p G 

The problems on which I shall essentially focus are the following ones: 

A. Find, if possible, optimal uniform decay estimates for the Fourier transform of 
the surface carried measure dp. 

B. If we denote by A t the averaging operator A t f(x) := j s f(x — ty)dp(y), determine 
for which exponents p the associated maximal operator 

Mf(x) := sup|A t /(x)| 

is bounded on L P (M. 3 ). For instance, if S is a Euclidean sphere centered at the origin, 
then M. is the spherical maximal operator studied first by Stein [61]. 

C. Determine the range of exponents p for which a Fourier restriction estimate 

1/2 



\f(x)\ 2 dp(x)) <C 



LP 



holds true. I shall explain how these questions can be answered in an (almost) complete 
way in terms of Newton polyhedra associated to the given surface S. All these results 
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are based on joint work with Ikromov, and in parts also with Kempe [39], [40], [37], 
[38] [41], [42]. 

Problem A is a classical question about estimates for oscillatory integrals, and there 
exists a huge body of results on it, in particular for convex hypersurfaces (some refer- 
ences, also to higher dimensional results, will follow). 

The other two problems had first been formulated by Stein: the study of maximal av- 
erages along hypersurfaces has been initiated in Stein's work on the spherical maximal 
function [61], and also the idea of Fourier restriction goes back to him (cf. §5). Their 
great importance for the study of partial differential equations became clear through 
Strichartz' article [63], and in the PDE-literature their dual versions are often called 
Strichartz estimates. 

Of course, these and related problems, such as the question of optimal sub-level 
estimates and integrability indices, make perfect sense also in higher dimensions (and 
possibly higher co-dimension), but only partial answers are known to most of these 
problems in this case. The reason for this is that strong information on the resolution 
of singularities is available for analytic functions of two variables, for instance by means 
of Puiseux series expansions of roots, whereas the situation for multivariate functions of 
more than two variables is substantially more complex. Nevertheless, there has been a 
lot of progress on the question as to how to construct more elementary and " concrete" 
resolutions of singularities for real analytic multivariate functions, giving more detailed 
information than what Hironaka's celebrated theorem [35] on the resolution of singular- 
ities would yield, for instance in work by Bierstone and P.D. Milman [7], [8], Sussman 
[64], Parusihski [53], [54], Greenblatt [28], [30], Collins, Greenleaf and Pramanik [15], 
among others. These techniques have already led to a very good understanding of, 
for instance, the sublevel estimation problem for slices of the surfaces in direction to 
the Gaussian normal, and the related determination of critical integrability indices, in 
independent work by Greenblatt [28], [30], and also Collins, Greenleaf and Pramanik 
[15], by rather different methods. Yet another approach has been developed by Denef, 
Nicaise and Sargos [21] in order to estimate oscillatory integrals in higher dimensions. 

These encouraging results give rise to the hope that eventually, substantial progress 
might also be possible in higher dimensions for the deeper problems A - C, but at 
present, only partial answers are available, for restricted classes of surfaces (see, e.g., 
[31], [33]). What makes these problems indeed a lot harder than the sub- level set 
problem is that it will surely require estimates of oscillatory integrals with phase func- 
tions depending also on small parameters, not just a fixed phase function. This leads 
to the fundamental problem of stability of the estimates in Problem A under small 
perturbation; I shall briefly come back to this later. 

Before returning to the two-dimensional case, let me address some results on the 
three problems A - C which have been obtained for particular classes of hypersurfaces 
also in higher dimensions. 
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The case of convex hypersurface of finite line type has been studied quite intensively, 
for instance in work by Randol [57], Svensson [65], Schulz [59], Bruna, Nagel and 
Wainger [12] and Cowling, Disney, Mauceri and myself [17] concerning problem A, 
and for instance by Nagel, Seeger and Wainger [52] and Iosevich, Sawyer and Seeger 
[46] on problem B. Various particular classes of non-convex hypersurfaces have been 
examined too, for instance in work by Iosevich and Sawyer [44], Cowling and Mauceri 
[18], [19], and by Sogge and Stein [60], where the method of damping (with powers 
of the Gaussian curvature) had been introduced in order to derive partial answers to 
problem B for very general hypersurfaces in R™. 

However, from now on I shall concentrate on the two-dimensional case, and give only 
occasionally a few references to work in higher dimensions. Since all three problems 
can be seen as "classical" by now, there is an abundance of literature associated to 
them, so that it would seem impossible to give due credit to everyone who has made 
contributions, and I apologize in advance to everyone whose work is not mentioned, 
due to lack of space or my personal ignorance. Many further references can be found 
in the cited articles, and for a more detailed account of the state of the theory and 
its historical development until roughly a decade ago, I refer, for instance, to Stein's 
monograph [62]. 

It is obvious that all three problems A - C can be localized to sufficiently small 
neighborhoods of given points x° on S. Observe also that the problems A and B are 
invariant under translations and rotations of the ambient space, so that we may replace 
the surface S by any suitable image under a Euclidean motion of R 3 . 

We may thus assume that x° = (0, 0, 0), and that S is the graph 



of a smooth function defined on a sufficiently small neighborhood f2 of the origin, 
such that 



The situation is different for problem B, since dilations do not commute with trans- 
lations, so that we are only allowed to work with linear transformations of the ambient 
space. We shall therefore study the maximal operator M. under the following transver- 
sality assumption on 5*. 

Assumption 1.1. The affine tangent plane x + T X S to S through x does not pass 
through the origin in R 3 for every x G S. Equivalently, x ^ T X S for every x G S, so 
that ^ S, and x is transversal to S for every point x G S. 

Notice that this assumption allows us to find a linear change of coordinates in R 3 
so that in the new coordinates S can locally be represented as the shifted graph of a 
function as before, more precisely, 



(1.1) 



S = {(x 1 ,X2,(f>(x 1 ,x 2 )) : (x 1 ,x 2 ) G fl} 



(1.2) 



0(0,0) = 0, V0(O,O) = 0. 



(1.3) 



S = {(xi, x 2 , 1 + (j){x u x 2 )) : (x ± , x 2 ) G Q}, 
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where <p satisfies again (1.2). 

Our transversality assumption is natural in this context. Indeed, various examples 
show that if it is not satisfied, then the behavior of the corresponding maximal function 
may change drastically. 

Observe also that if <fi is flat, i.e., if all derivatives of <fi vanish at the origin, and 
if p(x°) > 0, then it is well-known and easy to see that the maximal operator Ai is 
L p -bounded if and only if p = oo, so that this case is of no interest. In a similar way, 
also the problems A and C will become of quite a different nature. We shall therefore 
assume that is non-flat, i.e., of finite type. Correspondingly, we shall always assume 
that the hypersurface S is of finite type, in the sense that every tangent plane has finite 
order of contact. Recall that in the study a convex hypersurface, a standard assumption 
used by many authors is that the surface is of finite line type, which means that every 
tangent line has finite order of contact, which is stronger than our assumption. 

The article will be organized as follows: 

In the next paragraph, I shall briefly review some basic concepts concerning New- 
ton polyhedra, adaptedness of coordinates and the notion of height, which have been 
introduced by Arnol'd (cf. [2], [3]) and his school, most notably Varchenko [66]. These 
concepts had originally been studied for real analytic functions 0, but as shown in 
[39], can be extended to smooth, finite type functions <p. Problems A to C will then be 
discussed subsequently in §3 — §5. 

Since many of our proofs are quite elaborate, my main goal will be to give a kind of 
guided tour through the basic structure of the proofs, hoping that this might also be 
helpful for a everyone interested in studying some of the proofs in more detail. 

2. Newton polyhedra, and adapted coordinates 

Let me first recall some basic notions from [66], [39]. If is given as before, consider 
the associated Taylor series 

oo 

(j)(x 1 ,X 2 )~ ^ C »i,a2 X T x 2 2 

01,02=0 

of centered at the origin. The set 

T(0) := {(a x ,a 2 ) G N 2 : c Ql , Q2 = — L_^^0(O,O) ^ 0} 

ai'.a2'- 

will be called the Taylor support of at (0, 0). We shall always assume that 

T(0) + 0, 

i.e., that the function is of finite type at the origin. The Newton polyhedron M(4>) 
of at the origin is defined to be the convex hull of the union of all the quadrants 
(cm, 012) + in M 2 , with (011,012) G T(<j>). The associated Newton diagram J\fd(<f>) in 
the sense of Varchenko [66] is the union of all compact faces of the Newton polyhedron; 
here, by a face, we shall mean an edge or a vertex. 
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Figure 1 . Newton polyhedron 



We shall use coordinates (^,£2) for points in the plane containing the Newton poly- 
hedron, in order to distinguish this plane from the (£1,0:2) - plane. 

The Newton distance, or shorter distance d = d(4>) between the Newton polyhedron 
and the origin in the sense of Varchenko is given by the coordinate d of the point (d, d) 
at which the bi-sectrix t\ = t 2 intersects the boundary of the Newton polyhedron. 

The principal face vr(0) of the Newton polyhedron of <p is the face of minimal dimen- 
sion containing the point (d,d). Deviating from the notation in [66], we shall call the 
series 



the principal part of <fi. In case that vr(0) is compact, pr is a mixed homogeneous 
polynomial; otherwise, we shall consider pr as a formal power series. 

Note that the distance between the Newton polyhedron and the origin depends on 
the chosen local coordinate system in which is expressed. By a local coordinate system 
at the origin I shall mean a smooth coordinate system defined near the origin which 
preserves 0. The height of the smooth function is defined by 



where the supremum is taken over all local coordinate systems y = (2/1,2/2) at the 
origin, and where d y is the distance between the Newton polyhedron and the origin in 
the coordinates y. 

A given coordinate system x is said to be adapted to if h(<f>) = d x . 

In [39] we proved that one can always find an adapted local coordinate system in two 
dimensions, thus generalizing the fundamental work by Varchenko [66] who worked in 
the setting of real-analytic functions <p (see also [56]). 





(ai,a2)G7r(</>) 



h(4>) := sup{rf y }, 
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Recall also that if the principal face of the Newton polyhedron N{4>) is a compact 
edge, then it lies on a unique "principal line" 

L := {(ti,t 2 ) G M 2 : Mi + k 2 *2 = 1}, 

with K±, n 2 > 0. By permuting the coordinates X\ and x 2 , if necessary, we shall always 
assume that K\ < k 2 . The weight k = (k±, k 2 ) will be called the principal weight 
associated to 0. It induces dilations 5 r (xx,X2) := (r Kl £i, r K2 x 2 ), r > 0, on M 2 , so that 
the principal part pr of is k- homogeneous of degree one with respect to these 
dilations, i.e., <f> pr (S r (xi,x 2 )) = r<p pT (xi,x 2 ) for every r > 0, and 

(2.1) d= = -*-. 

Ki + K 2 \K\ 

It can then easily be shown (cf. Proposition 2.2 in [39]) that pr can be factorized as 

M 

(2.2) <Mxi,x 2 ) = cx^a:^ \[{x q 2 - X^r, 

1=1 

with M > 1, distinct non-trivial "roots" A; G C \ {0} of multiplicities 6 N\ {0}, and 
trivial roots of multiplicities z^i, v 2 G N at the coordinate axes. Here, p and q have no 
common divisor, and k 2 / Ki = p/q. 

More generally, if k = (ki,k 2 ) is any weight with < K\ < k 2 such that the line 
L K := {(ti, ^2) G : + K 2 t 2 = 1} is a supporting line to the Newton polyhedron 
jV(0) of 0, then the n-principal part of 

(ai,a2)Gl/«; 

is a non-trivial polynomial which is K-homogeneous of degree 1 with respect to the 
dilations associated to this weight as before, and which can be factorized in a similar 
way as in (2.2). By definition, we then have 

(2.3) <j)(xi,x 2 ) = 4> K (xi,x 2 ) + terms of higher ^-degree. 

Adaptedness of a given coordinate system can be verified by means of the following 
criterion (see [39]): Denote by 

m(0 pr ) := ord 5 i0 pr 

the maximal order of vanishing of pr along the unit circle S" 1 centered at the origin. 
The homogeneous distance of a K-homogeneous polynomial P (such as P = pr ) is 
given by dh(P) ■= + ^2) — l/M- Notice that (dh(P) , dh(P)) is just the point of 

intersection of the line given by K±ti + n 2 t 2 = 1 with the bi-sectrix t\ = t 2 . The height 
of P can then be computed by means of the formula 

(2.4) h(P)=max{m(P),d h (P)}. 

In [39] (Corollary 4.3 and Corollary 2.3), we proved the following characterization of 
adaptedness of a given coordinate system: 



8 



D. MULLER 



Proposition 2.1. The coordinates x are adapted to <p if and only if one of the following 
conditions is satisfied: 

(a) The principal face tt ((f)) of the Newton polyhedron is a compact edge, andm(cj) pr ) < 

d{4>). 

(b) 7r(0) is a vertex. 

(c) 7r(0) is an unbounded edge. 

These conditions had already been introduced by Varchenko, who has shown that 
they are sufficient for adaptedness when <fi is analytic. 

We also note that in case (a) we have h(4>) = h((fi pT ) = d h ((p vr ). Moreover, it can be 
shown that (a) applies whenever 7r(0) is a compact edge and k 2 J K\ 4- in this case 
we even have m(0 pr ) < d(<f>) (cf. [39], Corollary 2.3). 

2.1. Construction of adapted coordinates. In the case where the coordinates 
(21,22) are not adapted to <p, the previous results show that the principal face 7r(0) 
must be a compact edge, that m := k 2 /ki G N, and that m(0 pr ) > d((f>). One easily 
verifies that this implies that p = m, q = 1 in (2.2), and that there is at least one, 
non-trivial real root 22 = \ix± of </> pr of multiplicity ni = m(0 pr ) > d(<f). Indeed, one 
can show that this root is unique. Putting b\ := A;, we shall denote the corresponding 
root 22 = biXi of pr as its principal root. 
Changing coordinates 

Vi ■= xi, y 2 := 2 2 - hx™, 
we arrive at a "better" coordinate system y = (7/1,7/2)- Indeed, this change of coordi- 
nates will transform <p pr into a function pr , where the principal face of pr will be a 
horizontal half-line at level t 2 = m(0 pr ), so that d((j) pT ) > d(<f)), and correspondingly 
one finds that d(<j)) > d((j>), if (f> expresses is the coordinates y (cf. [39]). 

Somewhat oversimplifying, by iterating this procedure, we essentially arrive at Var- 
chenko's algorithm for the construction of an adapted coordinate system (cf. [39] for 
details). 

In conclusion, one can show that there exists a smooth real-valued function ip (which 
we may choose as the so-called principal root jet of 0) of the form 

(2.5) ^(21) = C2™ + 0(x™ +1 ) 

with c 7^ 0, defined on a neighborhood of the origin such that an adapted coordinate 
system (yi,y 2 ) for <p is given locally near the origin by means of the (in general non- 
linear) shear 

(2.6) 2/1 := 21, y 2 := x 2 - ^(^l)- 
In these adapted coordinates, <fi is given by 



(2.7) 



(t> a {y) := 0(2/1,2/2 + ^(3/1))- 
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Figure 2. <p(xi,x 2 ) :— (x 2 - xf) n + x[ (£ > mn) 

Example 2.2. 

(f)(xi, x 2 ) := (x 2 - x™) n + x{ . 

Assume that i > mn. Then the coordinates are not adapted. Indeed, 4> pr {xi,x 2 ) = 
(x 2 — x™) n , d(4>) = l/(l/n + l/ (mn)) = mnjim + 1) and m(0 pr ) = n > d(<p). Adapted 
coordinates are given by y\ := £1,2/2 := x 2 — x™, in which is expressed by <fi a (y) = 

y 2 +vl 

Remark 2.3. An alternative proof of Varchenko's theorem on the existence of adapted 
coordinates for analytic functions <fi of two variables has been given by Phong, J. Sturm 
and Stein in [56] , by means of Pusieux series expansions of the roots of (f). 

Let us finally observe that when m = K 2 j K\ = 1 in the first step of Varchenko's 
algorithm, then a linear change of coordinates of the form yi = xi, y 2 = x 2 — biXi will 
transform <fi into a function </>. Since all of our problems A - C are invariant under such 
linear changes of coordinates, by replacing our original coordinates (x\,x 2 ) by (yi,y 2 ) 
and 4> by 0, we may in the sequel always assume without loss of generality that either 
our coordinates (xi,x 2 ) are adapted, or they are not adapted and 

(2.8) m = K 2 j ' K\ is an integer > 2. 

A linear, non-adapted coordinate system for which (2.8) holds true will be called linearly 
adapted to 0. 

3. Problem A: Decay of the Fourier transform of the surface carried 

MEASURE dfl. 

Observe first that in view of (1.1), we may write /i(£) as an oscillatory integral 

£(f ) = : J(f ) = f e -i(^i^)HixiH2x 2 ) v ^ dx ^ ^ e r3) 
Jn 
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where rj G C^°(Q). Since V0(O, 0) = 0, the complete phase in this oscillatory integral 
will have no critical point on the support of r\ unless |£i| + |£ 2 | *C [£3!, provided fl 
is chosen sufficiently small. Integrations by parts then show that = 0(|^| _JV ) as 
|£| -> 00, for every JVeN, unless |£i| + |£ 2 | < 

We may thus focus on the latter case. In this case, by writing A = —£3 and £j = 
SjX, j = 1, 2, we are reduced to estimating two-dimensional oscillatory integrals of the 
form 

/(A; s) := J e lXWxi > X2)+SlXl+S2X2) r 1 (x 1 ,x 2 ) dx 1 dx 2 , 

where we may assume without loss of generality that A > 1, and that s = (si, s 2 ) G K 2 
is sufficiently small, provided that rj is supported in a sufficiently small neighborhood 
of the origin. The complete phase function is thus a small, linear perturbation of the 
function <fi. 

If s = 0, then the function /(A; 0) is given by an oscillatory integral of the form 
f e lX ^ x ^r](x) dx, and it is well-known ([6], [4]) that for any analytic phase function 
defined on a neighborhood of the origin in W 1 such that 0(0) = 0, such an integral 
admits an asymptotic expansion as A — > 00 of the form 

00 n—l 

(3.1) ^^ %ifc (0)A-^log(Ay, 

k=0 j=0 

provided the support of rj is sufficiently small. Here, the rj, form an increasing sequence 
of rational numbers consisting of a finite number of arithmetic progressions, which 
depends only on the zero set of 0, and the aj^ are distributions with respect to the 
cut-off function rj. The proof is based on Hironaka's theorem. 

Let us come back to the case n = 2. Following [66] (with as slight modification), we 
next define what we like to call Varchenko's exponent v{<fi) G {0, 1} : 

If there exists an adapted local coordinate system y near the origin such that the 
principal face vr(0 a ) of 0, when expressed by the function a in the new coordinates, is a 
vertex, and if h{4>) > 2, then we put v{4>) ■= 1; otherwise, we put v(4>) '■= 0. We remark 
[40] that the first condition is equivalent to the following one: If y is any adapted local 
coordinate system at the origin, then either ir(<f) a ) is a vertex, or a compact edge and 

m(r p r) = d((p a )- 

Varchenko [66] has shown that the leading exponent in (3.1) is given by r = l/h(4>), 
and is the maximal j for which aj t k{4>) 7^ 0. This implies in particular that 

(3.2) |/(A; 0)| < CA~mW log(Af w , A > 1, 

and this estimate is sharp in the exponents. Subsequently, Karpushkin [47] proved that 
this estimate is stable under sufficiently small analytic perturbations of (analogous 
results are known to be wrong in higher dimensions [66]). In particular, we find that 
J(A; s) satisfies the same estimate (3.2) for \s\ sufficiently small, so that we obtain the 
following uniform estimate for fi, 

(3.3) \m\ <c(i + ieir^ io g (2 + i£ir w , eet 3 , 
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provided the support of p is sufficiently small. 

In [40], we proved, by a quite different method, that Karpushkin's result remains 
valid for smooth, finite type functions <j), at least for linear perturbations, which led to 
the following 

Theorem 3.1. Let S = graph(0) be as before, and assume that (p is smooth and of 
finite type. Then estimate (3.3) holds true provided the support of p is sufficiently 
small. 

The special case where £ = (0, 0, £3) is normal to S at the origin is due to Greenblatt 
[29]. 

One can also show that this estimate is sharp in the exponents even when is not 
analytic, except for the case where the principal face n((p a ) is an unbounded edge. 
Indeed, if 7r(0 a ) is compact, then [40] 

/(A; 0) x CA"mW log(A)^ 

as A — > +00, where C is a non-zero constant. However, when 7r(0 a ) is unbounded then 
the following examples, due to Iosevich and Sawyer [45], show that there may be a 
different behavior in general: Let (p(xi,x 2 ) : = x\ + e -1 /' 211 '", with a > 0; then 

|/(A;0)| - A1/2lQ 1 gA1/Q as A -)• +00, 

whereas v{4>) = 0. 

It appears likely that even the full analog of Karpushkin's theorem, i.e., stability of 
estimate (3.2) under arbitrary small, smooth perturbations, holds true in two dimen- 
sions. 



3.1. Outline of some main ideas of the proof. Our proof of Theorem 3.1 is based 
on a certain decomposition scheme for the given surface S, related to the Newton 
polyhedra of <fi respectively a , which has been inspired by the work of Phong and Stein, 
in combination with re-scaling arguments. Since the same type of decompositions plays 
an important role also in the more involved proofs of our results on the other problems 
B and C (which require also further, more subtle refinements of them), I shall outline 
some of the major ideas used in the proof subsequently. 

We first notice that the case where h(<p) < 2 is covered by Duistermaat's work [24] 
(notice that Duistermaat proves estimates of the form (3.3) without the presence of a 
logarithmic factor log(2+ |£|), even for a wider class of phase functions). Note that the 
required estimates can also be derived easily from the normal forms in Theorem 5.8. 

I shall therefore subsequently assume that 

h := h(4>) > 2. 
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In many situations, one can reduce the problem to a one-dimensional one by means 
of van der Corput's lemma, respectively the following (not quite straight-forward) con- 
sequence of it, whose formulation goes back to J. E. Bjork (see [22]) and G. I. Arhipov 

[!]■ 

Lemma 3.2. Assume that f is a smooth real valued function defined on an interval 
I C M. which is of polynomial type n > 2 (n G N), i.e., there are positive constants 
Ci, c 2 > such that 

Ci < |/^( s )| < c 2 f or every s G /. 
i=2 

Then for A G E, 

r e iX ^g(s)ds\<C\\g\\ CH i ) (l + \X\)- 1/n , 
where the constant C depends only on the constants c\ and c 2 . 

3.1.1. The case where the coordinates are adapted to <p. Let us write d = d(<f)), and 
recall that here h = d, since the coordinates are adapted. Assume for instance that 
the principal face 7r(0) is a compact edge. By decomposing M 2 into the half-spaces 
M."± := M x M.±, we may also assume that the integration in J(£) takes place over one of 
these half-spaces only, say, IR+. Let k be the principal weight, with associated dilations 
5 r (xi,x 2 ) = (r Kl xi,r K2 X2). Recall also that then pr = (p K is 5 r -homogeneous of degree 
1. 

We fix a suitable smooth cut-off function \ on E 2 supported in an annulus A on 
which |x| ~ 1, such that the functions Xk '■= X $2 k form a partition of unity, and then 
decompose 

oo 
k=ko 

where 

J k (£):= [ e- l(?30(x)+ « lXl+6x2) r/(x)xfc(a;)dx. 



Scaling by 5 2 -k, we see that 

(3.4) J fc (0 = 2 / e V J V (5 2 - k (x))x(x)dx, 

Jr 2 + 

with fc (x) := 2 k (f)(5 2 -kx). Notice that in view of (2.3), 

4> k (x) = (f) K (x) + error term. 

We claim that given any point x° G A, we can find a unit vector e G M 2 and some 
j G N with 2 < j < h such that d J e (j) K (x ) ^ 0, where d e denotes the partial derivative 
in direction of e. 
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Indeed, if V<j) K (x ) 7^ 0, then the homogeneity of <j) K and Euler's homogeneity relation 
imply that rank (D 2 <j) K (x )) > 1. Therefore, we can find a unit vector e G M 2 such that 
dg4> K {x°) ^ 0. And, if V0 K (x°) = 0, then by Euler's homogeneity relation we have 
4> K (x°) = as well. Thus the function <p K vanishes in a; of order j > 2. But, in view of 
Proposition 2.1, then j < m(<f) pi .) < d = h, which verifies the claim. 

For k > ko sufficiently large we can thus apply van der Corput's lemma to the 
integration along lines parallel to the direction e in the integral defining near the 
point x°. Applying Fubini's theorem and a partition of unity argument, we thus obtain 



However, since we are assuming that 7r(0) is a compact edge, we have that = 
d(<fr) = h, and moreover M < h. This implies \k\M < 1. Since we have seen that here 
v{4>) = 1 if and only if M = m(0 pr ) = h, i.e., if and only if M\k\ = 1, we obtain 
estimate (3.3). 

3.1.2. The case where the coordinates are not adapted to 0. Here, in a first step, we 
may reduce to a narrow neighborhood of the principal root. 

Indeed, away from the principal root of <p w , we can argue in the same way as before, 
since the multiplicity of any real root of pr different from the principal root is bounded 
by d < h. I.e., we can reduce to a narrow K-homogeneous neighborhood of the curve 
£2 = bi%Ti °f the form 



say by means of a function pi(x) := Xo(( x 2 — biX™) / (ex™)) , where Xo is a suitable 
smooth bump function supported in the interval [—1,1] and e > is sufficiently small. 
I.e., in place of /(£), it suffices to estimate J Pl (£), where we write 



if x is an Y integrable function. 

Second Step: Domain decomposition into "homogeneous" domains D\ and transition 
domains E\. In oder to study the contribution J P1 (£) by the domain (3.7), we change 
to the adapted coordinates y and essentially re-write 



14(01 < L7|M| C 3 (K2) 2- fe N (1 + 2- fc |£ 3 |)- 1/j 
(3.5) < C|M| C3(R2) 2-^1(1 + 2- fc |el)- 1/M , 




(3.7) 



x 2 — bix™\ < ex 



.in 




(3.8) 
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Figure 3. Edges and weights 

where fj and Xo have properties similar to 77, respectively Xo- 

Let us denote the vertices of the Newton polyhedron Af((f> a ) by (A[, B[), I — 0, ... ,n, 
where we assume that they are ordered so that Ai-% < Ai, I — 1, ... ,n, with associated 
compact edges given by the intervals 7/ := -Bz_i), (Ai, Bi)], I = l,...,n. The 

unbounded horizontal edge with left endpoint (A n ,B n ) will be denoted by 7 n +i- To 
each of these edges 7;, we associate the weight n l = (k[, k 1 2 ), so that 7; is contained in 
the line 

U ■= {{h,t 2 ) G M 2 : K[h + K l 2 t 2 = 1}. 
For I = n + 1, we have < +1 := 0, k% +1 = l/B n . We denote by 

a t := ^j, I = 1, . . . ,n 

the reciprocal of the slope of the line L[. For I — n + 1, we formally set a n+ i := 00. 

If I < n, the ^'-principal part of (fi a corresponding to the supporting line Li is of 
the form 

(3.9) rAv) = ciyf 1 -^ 1 n (v2 - c?y?) Na 

a 

(cf. [38]). 

Remark 3.3. When <fi is analytic, then this expression is linked to the Puiseux series 
expansion of roots of (f> as follows [55] (compare also [39]): 
We may then factorize 

<p a (yuy2) = U( yi ,y 3 )y?y? J] (y 2 - r( Vl )), 

r 

where the product is indexed by all non-trivial roots r = r(y 1 ) of <p a (which may also 
be empty) and where [7(0, 0) 7^ 0. Moreover, these roots can be expressed in a small 
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a i—<*p-i 



neighborhood of as Puiseux series 

Kvi) = <%v? + Ci?y? x + ■■■ + €Z p y* lv ' 4p + ■ ■ 

where 



a\—OL p -\p _j_ ai—ap_i7 



ai-ap_i ai-ap_a 

wit/i strictly positive exponents a^..j ap_1 > and non-zero complex coefficients c^.. ; " p 7^ 
0, and where we have kept enough terms to distinguish between all the non-identical 
roots of <f> a . The leading exponents in these series are the numbers 

a\ < a 2 < ■ ■ ■ < a n . 

One can therefore group the roots into the clusters of roots [I], I — 1, ... ,n, where the 
I 'th cluster [I] consistes of all roots with leading exponent a;. 
Correspondingly, we can decompose 

n 

<p a (yi,V2) = U(y u y 2 )yl^l[^ ll] (y 1 ,y 2 ), 

1=1 

where 



% (2/1,2/2) := JJ^ 2 ~ r (^))- 



re 



Observe the following: If S l s (xi,x 2 ) = (s Kl Xi, s K2 Xi), s > 0, denote the dilations 
associated to the weight n l , and if r 6 [li] is a root in the cluster [li], then one easily 
checks that for y = (2/1,2/2) i n a bounded set we have 5 l s y2 = s K ^y 2 and r(5 l s yi) = 
s a h K icf i 1 yl h (1 + 0(s e )) as s — > 0, for some e > 0. Consequently, 

r -s^c^ 1 , if/i</, 

tflfa - r(^) = (1 + 0(s £ )) ^ s 4 (Ka _ c ^), if /! = / , 

t s K 2?/2, if h > l 

This shows that the /^-principal part of <p a is given by 

(3.10) r Kl = Cty^^^y^^ ]J(y 2 - <f 

where A^ )Ctl denotes the number of roots in the cluster [I] with leading term c^y® 1 . A 
look at the Newton polyhedron reveals that the exponents of y\ and y 2 in (3.10) can 
be expressed in terms of the vertices (A,-, Bj) of the Newton polyhedron: 

ii<Z ii>( 
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Notice also that 

n 

ai 1=1 

Comparing this with (3.9), the close relation between the Newton polyhedron of a 
and the Pusieux series expansion of roots becomes evident, and accordingly we say 
that the edge ji := [(Ai_±, Bi_i), (Ai, B{)] is associated to the cluster of roots [I]. 

Next we choose the integer Iq > 1 such that 

ai < • • • < Qfo-i < m < ai < ■ • • < ai < ai + i < ■ ■ ■ < a n . 

Since the original coordinates x were assumed to be non-adapted, the vertex (Aj _i, Bi ^i) 
will lie strictly above the bisectrix, i.e., Ai -\ < B io _x, . 

Let us consider the case where the principal face of the Newton polyhedron of <p a is 
a compact edge (the other cases require modified arguments). We choose A > Iq so 
that the edge 7a = Bx-i), (A\,B\)] is the principal face 7r(0 a ) of the Newton 

polyhedron of <p a (cf. Figure 3, where A = 3.) In the case where m(0p r ) = d(<p a ), it can 
easily be seen that by running Varchenko's algorithm one step further, we can pass to 
new adapted coordinates in which the principal face is a vertex. 

We therefore may assume that m(0p r ) < d(<p a ), so that v(4>) = 0. 

We shall narrow down the domain (3.7) to a neighborhood D\ of the principal root 
jet of the form 

(3.11) |x 2 -V(xi)| <N x x1\ 

where A^a is a constant to be chosen later. This domain is K A -homogeneous in the 
adapted coordinates y. 

To this end, following an idea from [55], we decompose the difference set of the 
domains (3.7) and (3.11) into the domains 

Di := {(x 1 ,x 2 ) : < \x 2 - ^(a?i)| < I = l , ■ ■ ■ , A - 1, 

and the intermediate domains 

Ei := {(x 1 ,x 2 ) : N l+1 xl l+1 < \x 2 - rj>(xi)\ < e^}, I = l , . . . , A - 1, 

as well as Ei _i := {(xi,x 2 ) : N^x"^ < \x 2 — ip(x{)\ < s^x™}. Here, the E\ > are 
small and the iVj > are large parameters to be chosen suitably. 
Observe that the domain 

Df := {{y u y 2 ) : e^ 1 < \y 2 \ < N iy ?} 

corresponding to Di is ^'-homogeneous in the adapted coordinates y given by (2.2), and 
contains the cluster of roots [/] , if Si and iVj are chosen sufficiently small, respectively 
large, while the domain Ef corresponding to Ei can be viewed as a domain of transition 
between two different homogeneities. 
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The contributions by the domain Di to J Pl (£) can again be estimated in a similar 
way as we did in the adapted case, by using dyadic decompositions and subsequent 
re-scalings by means of the dilations S l r associated to the weight k 1 . 

More precisely, the corresponding fc'th term will be given by the following analogue 
of (3.4) 

fc | K i| f ^(2- fe 60 fc (y)+6^(2- fcK iyi)+2- feK i €l y 1 +2- fc 4 52 y 2 j 



^(0 = 2-*!" I / e V Jrj(5 t 2 . k y) X (y)dy, 
Jm 2 + 

where 4> k (y) = <fi a K i{y) + error term. 

Notice, however, that since 1 — mn\ > k 2 — mK x > 0, the contribution of the non- 
linearity ip to the complete phase of the corresponding oscillatory integrals may be 
large, compared to the term containing fc , so that we are only allowed to apply van 
der Corput's estimate to the integration with respect to the variable 2/2 if we want to 
reduce to one- dimensional oscillatory integrals! 

We therefore need a control on the multiplicities of roots of d 2 4> a i at points y° in the 
corresponding annulus A not lying on the y\ axis (note that the latter corresponds to 
the principal root jet in the coordinates y). Indeed, it can be shown (cf. Propostion 2.3 
(b) in [38]) that these multiplicities are bounded by dh{<p a i) — 2, where dh{(t> a K i) denotes 
the homogeneous distance of (f> a u and it is evident from the geometry of the Newton 
polyhedron of 4> a that dh(<j> a i) < d(<p a ) = h, so that for every point y° in A H D\ there 
is some j G {2, . . . , h — 1} such that 



As for the contributions by the domains Ei, here we perform a separate dyadic 
decomposition in both variables y\ and 1/2, so that we geometrically decompose E\ into 
dyadic rectangles of size 2~- ? x 2 _fc , and then re-scale in both variables so that these 
rectangles become the standard cube, say, [1,2] x [1,2]. 

Now, the estimates of Section 11 in [38] show that the phase functions that one 
obtains after these re-scalings satisfy the estimate 

dl4>lk{y Q ) ^ for every y° G [1, 2] X [1, 2]. 

Since h > 2, this clearly suffices to obtain the necessary order of decay of the Fourier 
transform of these dyadic pieces. Moreover, scaling back to the original dyadic rect- 
angles, a careful analysis of the dependency of the corresponding estimates on the 
parameters j, k shows that it is indeed possible to sum theses estimates and obtain the 
same type of estimate for the contributions by the domains E\ as for the domains Di, 
even without logarithmic factor. 

Note that so far, we have always been able to reduce our estimations to van der 
Corput's lemma, respectively, Lemma 3.2, i.e., to one- dimensional oscillatory integrals. 
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3.1.3. Third Step: Study of the contribution by the homogenous domain D\ containing 
the principal root jet. What remains to be estimated is the contribution of the domain 

(3.11) to J(£). As we shall see, the study of this domain will in certain cases require 
the estimation of genuinely 2-dimensional oscillatory integrals, and a reduction to the 
one- dimensional case is no longer possible. 

We recall that according to our convention 

(3.12) m(4> a pi ) < d(^) = K 

so that u(4>) = 0. In the adapted coordinates y, the domain (3.11) is given by 1 2/2 1 < 
N\y® x , and we can cover it by a finite number of /t A -homogeneous subdomains of the 
form \y 2 — cy\ x \ < s^yl* , where c G [— A^, N\], and where, for a given c, we may choose 
60 > suitably small. 

Writing if)(xi) = x™u(xi), with a smooth function u satisfying oj(0) 7^ 0, we can 
thus reduce to estimating oscillatory integrals of the form 

(3.13) J c (0= / e*Mp(*^f-) V (y)dy, 
with a phase function 

F(y, := &<P a (y) + 62/1 + ^yTMvi) + 62/2 

depending on £ e R 3 . 

Arguing in a similar way as in the case of adapted coordinates, and recalling that 
0p r = , we may again perform a dyadic decomposition and re-scale by means of the 
dilations <5 A , in order to write 



k=ko 



where 



(3.14) J k (0 = 2~^ k [ ^^pi^^^Uv) X(V) dy, 

1 v £ oUi ' 



where 



F k (y, s) := 4> a pr ( Vl , y 2 ) + s lVl + S 2 y?u>{2 ^ k Vl ) + s 2 y 2 + error, 
where s := (si, s 2 , S 2 ) is given by 

Sl : = 2 ( 1 -^) fc ii, S2 := 2^-^^, S 2 := 2^~ m ^ k s 2 . 

^3 £3 

Note that 2 < m < a\ = k 2 /k± and k ^> 1, so that j^] ^> |s 2 |, and that 

~ 1 and \y 2 - cy" x \ < e 
for y in the support of the integrand of </&(£)• Recall also that we are assuming that 

lei ~ 161- 

One is thus led to the estimation of oscillatory integrals depending on certain pa- 
rameters (here si, s 2 , S 2 ) which may have various relative sizes. 
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The case where IS2I > M for some sufficiently large constant M ^> 1 can easily be 
treated by means of van der Corput's lemma applied to the y\- integration, so let us 
assume that \S 2 \ < M. Then \s 2 \ <C 1, provided we have chosen k sufficiently large. 

We may also easily reduce to the case where 

(3.15) <^>p r (l,c) = dirA^c) = for 1 <j < h, 

for otherwise an integration by parts in y% (if j = 1) or a simple application of Lemma 
3.2 yields a suitable estimate as before. 

The case where c > can easily be reduced to the case c = by performing another 
change of variables 2/2^2/2 + cy1 x in the integral defining Jfc(£). Indeed, one can show 
that our assumption (3.15) implies that necessarily a\ = fi^/^i £ N (cf. Corollary 3.2 
(iii) in [38]), and one checks that the new coordinates are again adapted to </>. 

So, let us assume that c = 0. Then necessarily 0" (1,0) 7^ 0, for otherwise 0p r would 
have a root of multiplicity at least h at (1,0), which would contradict (3.12). 

Assuming without loss of generality that 0" (1,0) = 1, we can then write (compare 
[38], Subsection 9.1) 

0pr (2/1,2/2) = 2/fQ(2/l,2/2)+2/l, 

where Q is a K A -homogeneous polynomial such that Q(l, 0) 7^ 0, and where B > h > 2. 

Recall that l^l < M, so that | S2 1 <C 1. Moreover, the case where |si| > N for some 
sufficiently large constant N ^> 1 can easily be dealt with by means of an integration 
by parts in y 1 , so let us assume that |si| < N. 

It will then suffice to show that, given any point (s°, S®) G [— M, M] x [-N, N] and 
any point y\ ~ 1, there exist a neighborhood U of (s°, S%), a neighborhood V of and 
some a > 1/h such that we have an estimate of the form 

n— k\n\ 

(3.16) |J fe (0| < 



for every (si,^) G U, provided the function x i n the definition of J&(£) is supported 
in V and Eq and k are chosen sufficiently small, respectively large. Summing over all 
k, this will clearly imply an estimate as in (3.3), even without logarithmic factor. 

To this end, first notice that for (s\, S2) G U and k sufficiently large, the function 
Fk(y, s) can be viewed as a small C°°- perturbation of the function 

F P r (y) ■= y^Qivi, 2/2) + s% + S 2 tu(0)y? + y\. 

Thus, if V-Fp r (y±, 0) 7^ 0, then we obtain (3.16), with a — 1, simply by integration by 
parts. 

Assume therefore that 0) is a critical point of F pr . Then y\ is a critical point of 
the polynomial function 

g{ yi ) := 4 + ^(0)^ + ^, 
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which comprises all terms of F pr depending on the variable y\ only. Note that the 
exponents satisfy 1 < m < n, since n = 1/k\ > fi^/^i > m - It is then easy to see 
that g" and g 1 " cannot vanish simultaneously at the given point y®, so that there are 
positive constants c\, c 2 > and a compact neighborhood V of (y®, 0) such that 



Cl 



< Ifl^G/OI ^ c 2 for ever y 2/1 e ^- 



This implies an analogous estimate for the partial derivatives d[Fk(yi, y 2 , s) of i^, 
uniformly for (si, S2) G f/ and the range of 2/2 that we consider, provided we choose U 
and Eq sufficiently small. Applying the van der Corput type estimate in Lemma 3.2, 
we thus obtain estimate (3.16), with a = 1/3, so that we are done provided h > 3. 
Notice also that if g"{y^) 7^ 0, then by the same type of argument we see that (3.16) 
holds true with a — 1/2 > 1/h. 

We may thus finally assume that 2 < h < 3, and that g'(yi) = g"{y1) = 0. In this 
case we have 

1 k a 

h < 3 and -4 > m > 2, 



so that 1/4 < 9/2. 

Note that B < \j n\ is a positive integer, and h < B < 9/2, so that either I? = 4 or 
-B = 3. We translate the critical point (yf, 0) of F pT to the origin by considering the 
function 

F« r (z) := F pr (y? + z 1? z 2 ) - g{yl) = z*Q{y\ + z u z 2 ) + ^ (3) (y?) 4 + ---- 

It is easy to see that this function has height W 1 := h(F$ T ) given by h$ = 1 / 3 +i/b < ^ - 
We can therefore again apply Duistermaat's results in [24] to the corresponding part 
of the oscillatory integral and obtain estimate (3.16), with a = > l//i. Note 
here that the estimates in [24] are stable under small perturbations. 

4. Problem B: L p - boundedness of the maximal operator M associated 

TO THE HYPERSURFACE S. 

We next turn to Problem B. As has already been mentioned in the introduction, 
the first, fundamental result on this problem is due to Stein [61], who proved that for 
n > 3, the spherical maximal function is bounded on L p (R n ) for every p > n/(n — 1). 
The analogous estimate in dimension n = 2 turned to require even deeper methods 
and was later established by Bourgain [9]. These results became the starting point 
for intensive studies of various classes of maximal operators associated to subvarieties. 
Stein's monograph [62] is an excellent reference to many of these developments. From 
these early works, the influence of geometric properties of S on the validity of L p - 
estimates for the maximal operator M. became evident. For instance, Greenleaf [34] 
proved that Ai is bounded on L p (R n ) if n > 3 and p > n/{n — 1), provided S has 
everywhere non- vanishing Gaussian curvature. In contrast, the case where the Gaussian 
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curvature vanishes at some point is still wide open, and complete answers are available 
at present only for finite type curves in the plane [43], and finite type hypersurfaces in 
M 3 and p > 2 (cf. Theorem 4.1). 

Let me come back to our hypersurface S in IR 3 . Recall from the introduction that 
we assume that S satisfies the transversality condition of §1, so that, by localizing to 
a small neighborhood of a given point x° G S and applying a suitable linear change 
of coordinates, we may assume that x° = and that S is given as the graph S = 
graph(l + 0), where is a smooth, finite type function defined on a neighborhood of 
the origin and satisfying 0(0, 0) = 0, V0(O, 0) = 0. We then define the height of S at 
x° by 



It easily seen that this notion is invariant under affine linear changes of coordinates in 
the ambient space M 3 . Recall also that JAf(x) := sup t>0 \A t f(x)\, where A t denotes 
the averaging operator over the t-dilate of S given by 



We can then state our main result from [38] (Theorems 1.2, 1.3), which gives an almost 
complete answer to the question of L p -boundedness of the maximal operator M. when 
p > 2 : 

Theorem 4.1. Assume the hypersurface S satisfies the transversality Assumption 1.1. 

(i) If the measure pda is supported in a sufficiently small neighborhood ofx°, then 
M. is bounded on L P {E?) whenever p > max{/i(x°, S), 2}. 

(ii) If M. is bounded on L P {E?) for some p > 1, and if p(x°) > 0, then p > h(x°, S). 
Moreover, if S is analytic at x°, then p > h(x°, S). 

4.1. Related quantities: contact index and sublevel growth. In [44], Iosevich 
and Sawyer had discovered a very interesting connection between the behavior of maxi- 
mal functions Ai associated to hypersurfaces S (in IR n ) and an integrability index asso- 
ciated to the hypersurface. In order to describe this, assume that dn(x) := dist (H,x) 
denotes the distance from a point x on S to a given hyperplane H. In particular, if 
x° G S, then c? T a .o(x) := dist (x° + T x oS,x) will denote the distance from x G S to the 
affine tangent plane to S at the point x°. The following result has been proved in [44] 
in arbitrary dimensions n > 2 and without requiring Assumption 1.1. 

Theorem 4.2 (Iosevich-Sawyer). If the maximal operator M. is bounded on L P (R"), 
where p > 1, then 



h{x°,S) := h{(f>). 




t > 0. 



(4.1) 




for every affine hyperplane H in M. n which does not pass through the origin. 
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Moreover, it was conjectured in [44] that for p > 2 the condition (4.1) is indeed 
necessary and sufficient for the boundedness of the maximal operator M. on L p , at 
least if for instance S is compact and p > 0. 

Remark 4.3. Notice that condition (4.1) is easily seen to be true for every affine 
hyperplane H which is nowhere tangential to S, so that it is in fact a condition on 
affine tangent hyperplanes to S only. Moreover, if Assumption 1.1 is satisfied, then 
there are no affine tangent hyperplanes which pass through the origin, so that in this 
case it is a condition on all affine tangent hyperplanes. 

Moreover, it is not very hard to prove (cf. [38]) that if S is a smooth hypersurface 
of finite type in R 3 , then, for every p < h(x°, S), 



and every neighborhood U of x°. And, if S is analytic near x°, then (4.2) holds true 
also for p = h(x°, S). 

Notice that this result does not require Assumption 1.1. As an immediate conse- 
quence of Theorem 4.1, Theorem 4.2 and (4.2) we obtain 

Corollary 4.4. Assume that S C M 3 is of finite type and satisfies Assumption 1.1, 
and let x° G S be a fixed point. Moreover, let p > 2. 

Then, if S is analytic near x°, there exists a neighborhood U C S of the point x° 
such that for any p G C£°(U) with p(x°) > the associated maximal operator M. is 
bounded on L P (IR 3 ) if and only if condition (4.1) holds for every affine hyperplane H 
in M 3 which does not pass through the origin. 

If S is only assumed to be smooth near x°, then the same conclusion holds true, with 
the possible exception of the exponent p = h(x°, S). 

This confirms the conjecture by Iosevich and Sawyer in our setting for analytic S, 
and for smooth, finite- type S with the possible exception of the exponent p = h(x , S). 
For the critical exponent p = h(x°, S), if S is not analytic near x°, examples show that 
unlike in the analytic case it may happen that M. is bounded on L h<yX (see, e.g., 
[45]), and the conjecture remains open for this value of p. 

In view of these results, it is natural to define the uniform contact index 7 lt (x°, S) of 
the hypersurface S at the point x° G S as the supremum over the set of all 7 for which 
there exists an open neighborhood U of x° in S such that the estimate 



holds true for every affine hyperplane H in IR n . If we restrict ourselves in this definition 
to the affine tangent hyperplane H = x° + T x oS at the point x°, we shall call the 
corresponding index the contact index j(x°, S) of the hypersurface S at the point 
x° G S. Here, we shall always assume that p(x°) 7^ 0. Note that if we change coordinates 
so that x° = and S is the graph of <j) near the origin, where <fi satisfies (1.2), then the 



(4.2) 




00 for every p < h(x , S) 
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contact index is just the supremum over all 7 such that there is some neighborhood U 
of the origin so that |</>|~ 7 G ^(U), so that the contact index agrees with the critical 
integrability index of (p as defined for instance in [15]. 

A closely related quantity is the sublevel growth rate cx(0), defined as the supremum 
over all a > such that there is some constant C a > so that 



Indeed, it can easily be shown by means of Tchebychev's inequality (cf. [15]) that 



In analogy with Arnol'd's notion of "singularity index" [2], let us finally introduce 
the uniform oscillation index (3 u (x°,S) of the hypersurface S at the point x° G S as 
the supremum over the set of all (3 such that 



for all p supported in a sufficiently small neighborhood of x . If we restrict directions £ 
to the normal to S at x°, then the corresponding decay rate will be called the oscillation 
index (3(x°, S) at x°. 

Combining our results with results from [56] (compare also [27]), we easily obtain 
the following result for smooth, finite type hypersurfaces S in IR 3 : 

Corollary 4.5. Let x° G S C M 3 be a given point so that p(x°) > 0. Then 



true. Indeed, Varchenko had already observed in [66] (Example 3) that for graphs 
S of functions <fi of three variables the oscillation index (3(0, S) may differ from 
l/h(0, S). Moreover, it is easy to give examples of functions of three variables 
where the contact index 7 = 7(0, S) is strictly smaller than the oscillation index 
(3 = (3(0, S). For instance, this applies to 0(xi,x 2 ,x 3 ) := x\ — {x\ + where 
the method of stationary phase shows that (3 = 3/2, whereas the factorization 
4>(xi, x 2 , x 3 ) = (x 3 — \J x\ + X2) (X3 + \/x\ + x%) shows that 7 = 1 (compare 
counterexample 8.1 in [15] j. 

However, the proof of Theorem 1.6 in [31] shows that the oscillation index (3 
can only be different from the contact index 7 when I/7 is an odd integer. 

To the best of my knowledge, it is not known whether the identities (3 u {x°, S) = 
f3(x°, S) and 7 U (£°, S) = j(x°, S) may persist in higher dimensions. Indeed, this 
question represents a special case of a question posed by Arnol'd, namely whether 
the oscillation index of a given function <fi is semicontinuos in the following 
sense: 

Suppose 4>(x, s) is a phase function which depends on x in a small neigh- 
borhood of the origin in W 1 and some small parameters s G M. k . Is there is a 
small neighborhood of (0, 0) G M. n x IR fc so that the oscillation index of <f>(x, s) 



{x G U : \<f>(x)\ < e} < C a e a for every s > 0. 



a(ct>)= 1 (x ,S), ifx 



and S = graph(^). 



\pd*(o\<Cp(i + \t\ye, 



Remarks 4.6. 
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at x° is greater or equal to that of (j)(x, 0) at the origin for every (x°, s°) in this 
neighborhood? That question had been answered to the negative in dimensions 
n > 3 by Varchenko [66], and to the positive for analytic functions depending 
on two variables by Karpushkin [47] . 

For linear perturbations of a given function <fr, our Theorem 3.1 shows that 
Arnol'd's conjecture holds true even for smooth, finite type functions <p, and 
it seems open whether this kind of stability may still hold true even in higher 
dimensions, since all counter examples to Arnol'd's conjecture which have been 
found hitherto have non-linear perturbation terms. 

Further important papers dealing with this stability question for general per- 
turbations are, e.g., [13], [56]. 

(b) Greenblatt [28], [30], and independently Collins, Greenleaf and Paramanik [15], 
have devised (quite different) algorithms of resolution of singularities which in 
principle allow to compute the contact index also in higher dimensions. 

(c) If p < 2, then examples (see, e.g., [46] j show that neither the notion of height 
nor that of contact index will determine the range of exponents p for which the 
maximal operator Ai is LP -bounded. Our work (in progress) on this case seems 
to indicate a certain conjecture how to express this range in terms of Newton 
polyhedra. Moreover, for certain surfaces this conjecture relates to fundamen- 
tal open problems in Fourier analysis, such as the conjectured reverse square 
function estimate for the cone multiplier (see, e.g., [49], [10],). 

4.2. A few hints on the proof of Theorem 4.1. Since the proof in [38] is quite 
involved, I shall here just try to indicate some of the main steps, grossly oversimplifying 
some of the arguments. 

In the first part of the proof, we basically follow the first two steps of the scheme 
that has been outlined in the preceding paragraph. Recall that in these steps, it had 
effectively been possible to reduce the problem to a one-dimensional one. The same 
applies basically to the estimation of the maximal operator A4, because in the first two 
steps, we have had a sufficiently good control on the multiplicity of roots by the height 
h. 

In these arguments, the following result plays an analogous role for problem B as 
the van der Corput type Lemma 3.2 played for problem A. 

Let U be an open neighborhood of the point x° G R 2 , and let pr G C°°(U, R) such 
that 

(4.3) d 2 m 0p r (x?,x°)^O, 
where m > 2. Let 

= 0pr + 0r, 

where <p r G C°°(U, R) is a sufficiently small perturbation. Denote by S £ the surface in 
R 3 given by S s := {(xi,x 2 , 1 + e<p(xi, x^)) '■ (^1,^2) G U}, with e > 0, and consider the 
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averaging operators 

A\f{x) := / f(x-ty)iP(y)da(y), 
Js £ 

where da denotes the surface measure and -0 G C£°(S e ) is a non- negative cut-off func- 
tion. Define the associated maximal operator by 

M e f{x) := sup 1^/(2)1- 

Proposition 4.7. Assume that pr satisfies (4.3) and that the neighborhood U of the 
point x° is sufficiently small. Then there exist numbers M G N, 5 > 0, such that 
for every <p r G C°°(U,M.) with ||<^v||c A/ < ^ an d every p > m there exists a positive 
constant C p such that for e > sufficiently small the maximal operator Ai s satisfies 
the following a priori estimate: 

(4-4) \\M £ f\\ p <C p e^\\f\\ p , feS(R 3 ). 

This result can be reduced to the one- dimensional case, i.e., the study of maximal 
functions associated to curves in the plane. Indeed, just consider the "fan" of all 
hyperplanes passing through the Xi-axis. This fan will fibre the given surface into a 
family of curves, and our dilations leave each of these planes invariant. 

The proof for the analogous result for plane curves then basically follows Iosevich's 
approach in [43] . The case where m = 2 is indeed the most difficult on; the cases where 
m > 3 can then easily be reduced to this case by means of dyadic decompositions and 
re-scalings. Note that if m — 2, then the corresponding maximal operator in the 
plane behaves in a similar way as the "circular" maximal function studied by Bourgain 
[9]. An alternative approach to Bourgain's one had later been given by Mockenhaupt, 
Seeger and Sogge [51], [50], based on suitable local smoothing estimates for classes of 
Fourier integral operators. This approach is stable under small perturbations of the 
given curve, which is exactly what is needed for our purposes. 

Indeed, in our applications of Proposition 4.7, we even need to replace (p(xi,X2) by 
(f)(xi,X2 — ^ e (xi)), where the function ip £ may blow up like 0(e~ 5 ) for some 5 G [0, 1[, 
and it turns out that this is still admissible, i.e., the estimate (4.4) remains valid also 
in this case. 

Now, if the coordinates are adapted to <j) (a similar argument will apply to the first 
step when the coordinates are not adapted), then, in analogy with the decomposition 
J(Q = J2kLk we can dyadically decompose the surface S as in Subsection 3.1.1, 

and accordingly decompose 

oo 

A t f(y,y 3 ) = J2^ k tf(y,y3), 

k=ko 

where one finds that 

A k J(y,y 3 ) = 2"^ / f(y - tS 2 - k (x), y 3 - t(l + 2- h cf> k (x)) v(S 2 -hx) X (x) dx, 
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with (f) k is as before. We denote by A4 k the maximal operator associated to the aver- 
aging operators A k . Observe next that the scaling operators T k , defined by 

T k f(y,y 3 ) := 2^/(^2/), 2/ 3 ), 
act isometrically on L P (IR 3 ), and 

(T- k A k T k )f(y,y 3 ) = 2- k W [ f (y - tx,y 3 - t(l + 2~ k ( j> k (x))) V (5 2 - k (x))x(x) dx, 

so that we are reduced to estimating the maximal operator associated to this family 
of averaging operators. This, in return, can be accomplished by means of Proposition 
4.7, where we choose e := 2~ k . The resulting estimates can then be summed over k, 
and we arrive at the desired estimation for Ai in this case. 

Assume next that the coordinates x are not adapted to <p. Eventually, we then again 
arrive at the situation studied in the third step (compare Subsection 3.1.3), i.e., we 
need to estimate the contribution of the domain D\ to the maximal operator Ai. Again, 
one can indeed reduce to a small subdomain of the form 

\x2 — ip(xi) — cx\ x \ < Sqx1 x , with xi > 0, 

and also assume that (3.15) holds true at the point (1, c) = (1, 0). Recall that we may 
then write 

0^(2/1,2/2) = 2/2^(2/1, 2/2) + 2/", 

where Q is a K A -homogeneous polynomial such that Q(1,0) 7^ 0, and B > h > 2. 
Again, in this situation, a reduction to the one-dimensional case is no longer possible. 

It turns out that this is the most difficult case, which requires various further ideas. 
The heart of the matter are in fact rather precise estimations of certain classes of 
two-dimensional oscillatory integrals depending on small parameters (cf. §5 in [38]). 

Indeed, in order to understand the behavior of <fi a as a function of y 2 , for y\ fixed, 
we decompose 

r(yi,y2) = r(yi,o) + e( yi ,y 2 ), 

and write the complete phase $ for J(£) in adapted coordinates (1/1,1/2) in the form 

i ? (y,O = (60°(?/i,o) + &J/1 + &K3/1)) + (&0(z/i,j/2) + 62/2). 

Notice that 

a (j/!,o) = jtfp(yi), V(y0 = vTMvi), 0AVUV2) = viQ(vi,V2), 

where 9 K x denotes the K A -principal part of 9, and where p and u do not vanish at 
2/i=0. 

Now, by means of a dyadic decomposition and re-scaling using the K A -dilations 
{$r}r>o we would like to reduce our considerations as before to the domain where 
yi ~ 1. In this domain, \y 2 \ 1, so that 9 K x{y) ~ y 2 Q(yi,0). What leads to problems 
is that the "error term" 9 r := 9 — 9 K \, which consists of terms of higher K A -degree than 
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6 K \, may nevertheless contain terms of lower 2/2-degree lj < B of the form Ojy^Vx ■> P ro ~ 
vided rij is sufficiently large (this corresponds to a "fine splitting of roots" of 9 when 
is analytic). After scaling the fc-th dyadic piece in our decomposition by 8^_ k in order 
to achieve that y\ ~ 1 and | g/ 2 1 < Eq, such terms will have small coefficients compared to 
the one of y 2 Q{yi, IJ2), but for 1 2/2 1 very small they may nevertheless become dominant 
and have to be taken into account. 

In order to resolve this problem, we apply a further domain decomposition by means 
of a suitable stopping time argument, into homogeneous domains D' e and transition 
domains E' £ , oriented, in some sense, at the level sets of d 2 <p a , which in return again are 
chopped up into dyadic respectively bi-dyadic pieces. After re-scaling, the contributions 
of these pieces to the maximal operators can eventually be estimated by means of 
oscillatory integral technics in two variables. 

More precisely, it turns out that what is needed are uniform estimates for various 
classes of oscillatory integrals of the form 

J(X,a,8):=[ e iXF(x ' a ' 5) a(x,5)dx, (A > 0), 



with a phase function F of the form 

F(xi,x 2 ,o-,5) := fx(xx,S) + crf 2 (x 1 ,x 2 ,8), 

and an amplitude a defined for x in some open neighborhood of the origin in IR 2 with 
compact support in x. The functions /1, f 2 are assumed to be real- valued and depend, 
like the function a, smoothly on x and on small real parameters 81, . . . ,8 V , which form 
the vector 5 := (8±, . . . , 8 V ) G W. cr denotes another small real parameter. 

With a slight abuse of language we shall say that ijj is compactly supported in some 
open set U C IR 2 if there is a compact subset K C U such that supp ■?/>(•, 8) C K for 
every 8. 

To give an idea as to which type of oscillatory integrals we need to estimate, let 
us remark that the most difficult instance are "oscillatory integrals of degenerate Airy 
type", for which we have the following result: 

Theorem 4.8. Assume that 

1^(0, 0)\ + \d 2 Ji (0,0)| + 1^/1 (0,0)1 ^ and «W 2 (0, 0, 0) ^ 0, 

and that there is some m > 2 such that 

4/2(0,0,0) =0 forl = l,...,m-l and 9^/2(0, 0,0) ^0. 

Then there exist a neighborhood U C M 2 of the origin and constants e, e' > such 
that for any amplitude a which is compactly supported in U the following estimate 

(4-5) \J(X,a,8)\< ° 

holds true uniformly for \a\ + \8\ < s', where l m := | and c m := 1 for m < 6, and 

lm := 2(2m-Z) a7ld Cm := 2 f° T 111 - 6 ' 
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Observe that the order of decay 0(A -1 / 2-e ) in (4.5) is just what we need in order to 
apply the usual method to control the maximal operator on L 2 by means of Sobolev's 
embedding theorem applied to the scaling variable t > 0. 

Remark 4.9. If the Fourier transform of the surface carried measure fi decays more 
slowly than Od^l" 1 / 2 " 5 ), then one cannot directly apply Sobolev's embedding theorem 
in order to control the maximal function on L 2 . This happens in many situations when 
the Gaussian curvature of the hypersurface S vanishes. One method to overcome this 
problem is to use suitable damping factors in the amplitude of the corresponding oscil- 
latory integrals, for instance powers of the Gaussian curvature, and combine this with 
a suitable complex interpolation argument. This technique had been introduced in [60], 
and also been used in our first approach our problems in [37]. 

However, the choice of a suitable damping factor can become quite tricky a task, and 
we believe that the techniques in [38], which avoid damping factors and rather rely on 
suitable decompositions of the given surface in combination with re-scaling arguments, 
are simpler and more straight-forward. 

For an approach based on damping methods, we refer to [32], where most of those 
cases are treated by damping techniques which can be reduced to a one- dimensional 
problem. 

5. Problem C: Fourier restriction to the hypersurface S. 

We finally turn to the last problem, namely the determination of the range of expo- 
nents p for which an L p -L 2 Fourier restriction estimate 

(5.1) ( I \f\x)\ 2 d^x)) 1 ' 2 < C\\f\\ LPm , f G <S(M 3 ), 

J s 

holds true. 

The idea of Fourier restriction goes back to Stein, and a first instance of this concept 
is the sharp L p -L q Fourier restriction estimate for the circle in the plane by Zygmund 
[67], who extended earlier work by Fefferman and Stein [25] (see also Hormander [36], 
[14] for estimates on more general oscillatory integral operators). For subvarieties 
of higher dimension, the first fundamental result was obtained (in various steps) for 
Euclidean spheres S*™ -1 by Stein and To mas, who eventually proved that an LP-L 2 
Fourier restriction estimate holds true for S' n_1 , n > 3, if and only if p' > 2(2/ (n — 1) + 1) 
(cf. [62] for the history of this result). 

Even though I shall not pursue that problem here, let me briefly remind that a 
more general and even substantially deeper problem is to determine the exact range of 
exponents p and q for which an L p -L q restriction estimate 

( / \f(x)\Ua(x)) 1/9 <C\\f\\ LP{Rn) 

holds true for spheres. It is conjectured that this is the case if and only if p' > 2n/(n— 1) 
and p' > q(2/(n — 1) + 1), and there has been a lot of very deep work on this problem 
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by many mathematicians, including Bourgain, Wolff, Vargas, Vega, Katz, Tao, Keel, 
Lee, and most recently Bourgain and Guth [11], which has led to important progress, 
but the problem is still open in dimensions n > 3 and represents one of the major 
challenges in Euclidean harmonic analysis, bearing various deep connections with other 
important open problems, such as the Bochner-Riesz conjecture, the Kakeya conjecture 
and Sogge's local smoothing conjecture for solutions to the wave equation. I refer to 
Stein's book [62] for more information on these topics and their history until 1993, and 
various related essays by T. Tao which can be found on his webpage. 

Coming back to the restriction estimate (5.1) for our hypersurface S in M 3 , we 
begin with the case where there exists a linear coordinate system which is adapted to 
the function 0. For this case, a complete answer had been given in [40] (for analytic 
hypersurfaces, partial results had been obtained before by Magyar [48]) : 

Theorem 5.1. Assume that, after applying a suitable linear change of coordinates, the 
coordinates (xi,x 2 ) are adapted to 0. We then define the critical exponent p c by 

(5.2) p' c := 2h(4>) + 2, 

where p' denotes the exponent conjugate to p, i.e., l/p+ 1/p' = 1. 

Then there exists a neighborhood U C S of the point x° such that for every non- 
negative density p e C^°(U) the Fourier restriction estimate (5.1) holds true for every 
p such that 

(5.3) l<P<Pc- 

Moreover, if p(x°) ^ 0, then the condition (5.3) on p is also necessary for the validity 
of (5.1). 

In many cases, this result is an immediate consequence of Theorem 3.1 and Green- 
leaf's classical restriction estimate in [34], which I shall state here for the special case 
of hypersurfaces only. 

Theorem 5.2 (Greenleaf). Assume that < l^l -1 ^. Then the restriction esti- 

mate (5.1) holds true for every p > 1 such that p' > 2(h + 1). 

Indeed, observe that a problem arises only when Varchenko's exponent z/(0) equals 
one in (3.3). In this direct application of Greenleaf 's result yields only the range 

1 < P < Pc- 

To capture also the endpoint p = p c , recall from Subsection 3.1.1 that we had 
effectively decomposed the measure \i into a dyadic sum /x = YlkLko^k, where Hk = 
(Xk ® l)/i, and where /tfc(£) = Jk(C) is given by (3.4). Moreover, we had estimated 
Jfc(£) in (3.5), which in particular yields the estimate 

< c2-*i*i(i +2- fe ie 3 |)- 1/ft . 

The measures \x\. are supported in dyadic annuli of "radius" 2~ k , which are images 
of an annulus of radius of size one under the dilations 5 2 -k associated to the principal 
weight k, so that we cannot directly apply Greenleaf 's restriction estimate. However, it 
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is important to notice that our estimates for do not carry a logarithmic factor yet 

(that only arose for /}(£) in certain cases through summation over the k), and a simple 
re-scaling argument can then be applied to derive the following uniform restriction 
estimate for the family of measures /i^ : 

(5-4) J \f(x)\ 2 d^ k (x)<C 2 \\f\\l. 

Fix a cut-off function \ £ C^°(1R 2 ) supported in an annulus centered at the origin 
such that x = 1 on the support of Xi an d define dyadic frequency decomposition 
operators A' k by 

A>J(x) :=x(5 2k x')f(x',x 3 ), 

where we have written (xi,x 2 ) = x' . Then J \f(x)\ 2 dbtk(x) = J \A' k f(x)\ 2 dfik(x), and 
setting p := p c , we see that (5.4) yields in fact 

\hx)\ 2 d^(x)<C 2 \\ATft 
for any k > k . In combination with Minkowski's inequality, this implies 

1/2 / \ V2 / \ 1/2 



|/(x)| dfi(x 



) = (E / \fc)\ 2 dVL h {x)\ < [J2\\ A 'kf\\l 



c (E (/ i A fc/wr^) 



k>ko / \k>ko 

\ 2/p N 



, p/2\ V? 




) J £C 


(x: iai/(x)i 2 ) ' 




\k>k Q J 



LP 



since p < 2. Estimate (5.1) for p = p c follows thus by means of Littlewood-Paley 
theory. 

From now on, we shall therefore always make the following 

Assumption 5.3. There is no linear coordinate system which is adapted to <ft. 

According to our discussion in Subsection 2.1, we may then also assume that the 
coordinates x are linearly adapted to <fr, and that there are adapted coordinates y of 
the form y 1 = x x , y 2 = x 2 — 4>(xi), where 

(5.5) ip^) = x?u{x x \ with w(0) ^ and m > 2. 

4> a (y) will again denote <fi when expressed in these adapted coordinates, and we shall 
use the notions introduced for the study of the Newton polyhedron of <fi a in Subsection 
3.1.2. 

Consider the line parallel to the bi-sectrix 

A (m) := {(t,t + m + l) : t £ R}. 
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Figure 4. r- height 



For any edge 7; C Li := {(ti,t 2 ) G K 2 : K,[ti + K l 2 t 2 = 1} define hi by 

AHnL^^-m.Hl)}, 

i.e., 

1 + mK\ — hi 2 



(5.6) 



and define the restriction height, or short, r-height, of by 



/i r (0) := max(d, 



max 

{/=!,..., n+l:a;>m} 



/*/)■ 



Remarks 5.4. (a) For L m place of L\ and k in place of k 1 , one has m = k 2 / 1 ki 
and d — 1/(ki + k 2 ), so that one gets d in place of hi in (5.6). 
(b) Since m < aj, we have hi < + k' 2 ), hence h r (4>) < h(<j)). 

It is easy to see by Remark 5.4 (a) that the r-height admits the following geometric 
interpretation: 

By following Varchenko's algorithm (cf. Subsection 8.2 of [38]), one realizes that the 
principal line L of H{4>) is a supporting line also for the Newton polyhedron of 4> a , 
which intersects M{<p a ) in a compact face, either in a single vertex, or a compact edge. 
I.e., the intersection contains at least one and at most two vertices of J\f((j) a ), and we 
choose (v4; _i, -Bi -i) as the one with smallest second coordinate. Then Z is the smallest 
index / such that 7/ has a slope smaller than the slope of L, i.e., a; _i < m < ai . We 
may thus consider the augmented Newton polyhedron J\f r (<j) a ) of <ft a , which is the convex 
hull of the union of J\f(<j) a ) with the half-line L + C L with right endpoint (Ai _i, Bi Q _i). 
Then h r ((p) + 1 is the second coordinate of the point at which the line A(' m ) intersects 
the boundary of Af r (4> a ). 

The main result from [41], [42] then reads as follows. 
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Theorem 5.5. Let <fi 7^ be real analytic, and assume that there is no linear coordinate 
system adapted to <fi. Then there exists a neighborhood U C S of x° = such that for 
every non-negative density p G C£°(U), the Fourier restriction estimate (5.1) holds 
true for every p > 1 such that p' > p' c := 2h r (<f)) + 2. 

Remarks 5.6. (a) An application of Greenleaf 's result would imply, at best, that 
the condition p' > 2h(<p) + 2 is sufficient for (5.1) to hold, which is a strictly 
stronger condition than p' > p' c . 

(b) A. Seeger recently informed me that in a preprint, which regretfully had re- 
mained unpublished, Schulz [58] had already observed this kind of phenomenon 
for particular examples of surfaces of revolution. 

(c) It can be shown that the number m is well-defined, i.e., it does not depend on 
the chosen linearly adapted coordinate system x. 

Example 5.7. 

(f>(xi, x 2 ) := (x 2 — x™) n , n, m > 2. 

The coordinates (xi,x 2 ) are not adapted. Adapted coordinates are y\ := x\,yi := 
X2 — x™, in which is given by 

Here 

1 1 

mn n 

d := d((p) = ■ = — < n, 

K\ + k 2 m + 1 

and 

' _ / 2d + 2, if n<m + l, 
^ c 1 2n, if n > m + 1 . 

On the other hand, h := h((f)) = n, so that 2h + 2 = 2n + 2 > p' c . 

An analogous theorem holds true even for smooth, finite type functions 0, under an 
additional Condition (R) which, roughly speaking, requires that whenever the Newton 
diagram suggests that a root with leading term given by the principal root jet ip{x\) 
should have multiplicity B, then indeed such a root of multiplicity B does exist (this 
is a condition on the behavior of flat terms). Condition (R) is always satisfied when 
is real- analytic. 

For example, Condition (R) would hold true for 

<f>g{xi, X 2 ) = (X 2 ~x\- f{x x )f, 

for every flat smooth function f{x\) (i.e., /^(0) = for every j £ N). On the other 
hand, (R) is not satisfied for 

4>b(xi,x 2 ) := (x 2 - x\) 2 + f(xi), 

unless / vanishes identically. 
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I also like to mention that there is a more invariant description of the notion of 
r- height (cf. Proposition 1.9 in [41]), somewhat in the spirit of Varchenko's definition 
of height, but I refrain from stating it here since this would require the introduction of 
further, somewhat technical notions. 

5.1. Necessity of the condition p' > 2h r (<f)) + 2. In order to better understand 
the meaning of the notion of r- height, let me present the proof of the necessity of the 
condition p' > 2h r (4>) + 2 for the validity of the Fourier restriction estimate (5.1) when 
p(x°) 0. The proof will be based on a modified Knapp-type argument. 

Let 7; be any edge of Af((fi a ) with a>i > m, and choose the weight k 1 such that 7; lies 
on the line Li given by K,\ti + K l 2 t 2 = 1. Consider the region 

D a £ := {yeR 2 : \ Vl \ < e**, \y 2 \ < e^}, e > 0, 
in adapted coordinates y. In the original coordinates x, it corresponds to 

D £ := {x eM 2 : \xi\ < e K \ \x 2 - ip(xi)\ < e K ' 2 }. 
Assume that e is sufficiently small. Since 

rie^e^yi) = e + 0{e 5 ) ^ 

for some 5 > 0, we have that \4> a (y)\ < Ce for every y G D", i.e., 
(5.7) 10(^)1 < Ce for every x e D £ . 

Moreover, for x e D e , 

\x 2 \ <£ 4 + |V(^i)| <£^+e mK K 
Since m < ai — K l 2 / k\, we find that 

\x 2 \<e m <, 

so that we may assume that D £ is contained in the box where \xi\ < e^' 1 , \x 2 \ < £ mre i. 
Choose f e such that 

fe(xi, x 2 , x 3 ) = xo (|r) Xo (-^r) Xo ) . 
Then by (5.7) we see that f E {xi, x 2 , 4>(xi, x 2 )) > 1 on D £ , hence, if p(0) ^ 0, then 

1/2 



\? £ \ 2 pda) > |£) 6 |V2 = e (4+4)/2 



Since ||/ e || p — e(( 1+m ) K i +1 )/p' ; we find that the restriction estimate can hold true only if 

p , S2 a+ r w+i =2ft[+2i 

Aj^ Ay 2 

where we recall that /ij = (1 + mK,\ — k\) / {k\ + k\)- 
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Notice that the argument still works if we replace the previous line Li by the line 
L associated to the weight k, and (p a K i by 0". Since here mK\ = k 2 , this leads to the 
condition p' > 2d + 2, so that altogether necessarily 

p' > 2max((f, max hi) + 2 = 2h r (<p) + 2. 

l:a l >m 

Q.E.D. 

5.2. Sufficiency of the condition p' > 2h r ((p)+2: I. Some key steps in the proof 
when h\ in ((p) > 2. Let us assume that we are working in linearly adapted coordinates 
x, so that d := d((f>) = h\ in := h\ in ((j)). 

In the preceding discussions of problems A and B, it had been natural to distinguish 
between the cases where h := h(<p) < 2 and where h > 2, since in the latter case, 
in many situations a reduction to a one-dimensional situation had been possible by 
means of the van der Corput type Lemma 3.2. For similar reasons, in the discussion of 
Problem C it appears natural to distinguish between the cases where d < 2 and where 
d > 2. In addition, when d > 2, it turns out that the case where d > 5 can be handled 
in a somewhat simplified way compared to the case where 2 < d < 5. 

We shall therefore assume in this section that d > 5. In a first step, we can again 
localize to the narrow /t-homogeneous subdomain (3.7) of the curve x 2 = b\x™ given 
by 

\x 2 - cx™\ < ex™, 

by means of a cut-off function p\. Indeed, the technique of proof that we used in the case 
of adapted coordinates can essentially be carried over to the domain complementary 
to (3.7) without major new ideas, since one can show that the Fourier transforms of 
the corresponding dyadic pieces pk of the measure p satisfy estimates of the form 

|/^(0| <C2^l(l + 2-^ 3 |)- 1/d . 
Notice also that h r := h r (<f)) > d. 

Let us again assume for instance that the principal face of the Newton polyhedron 
of cf) a is a compact edge. Using the same notation as in Section 3, we choose again 
A > /o so that the edge 7a = [(Aa-i, -Ba-i), {A\, B\)] is the principal face vr(0 a ) of the 
Newton polyhedron of <p a . 

In a second step, we again narrow down the domain (3.7) to the neighborhood D pr := 
D\ of the principal root jet given by (3.11), where 

\x 2 - ^(xi)| < N x x1 x , 

again by decomposing the difference set of the domains (3.7) and (3.11) into the do- 
mains 

Di := {(xi,x 2 ) : Six^ 1 < \x 2 - i>{xi)\ < Nixf}, I = l , . . . , A - 1, 
and the intermediate domains 

Ei := {(xi,x 2 ) : Ni +1 x" l+1 < \x 2 - i/j(xi)\ < £ix\ 1 }, I = l , . . . , A - 1, 
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as well as £^ _i := {(xi,X2) : Ni x 1 ° < \x% — if)(xi) \ < £ix™}. 

Contribution by the domains E\. Let us denote by fiE t the contribution of the 
transition domains E\ to the measure /x. We then decompose bi-dyadically w.r. to 
the adapted coordinates y as 

so that fij t k is supported where y\ — x± ~ 2~ J and 2/2 = ^2 — if>(xi) ~ 2~ fc . Observe that 
this a "curved rectangle" in our original coordinates x. In a similar way as in the case 
of adapted coordinates, we would like to localize to these curved rectangles by means 
of Littlewood-Paley theory in order to reduce to uniform restriction estimates for the 
family of measure /x,^, i.e., 

(5-8) / |/| 2 ^-, fc <C|U|| 2 

Js 

for p < p c . Clearly, because of the non-linearity ip(xi), this is not possible by means 
of Littlewood-Paley techniques in the variables X\ and x%, but it turns out the we can 
use the variables x\ and 23 to accomplish this. 
Indeed, one can show that 

a (l/) = c l V\ l 1h l + small error j on Ef, 

which in return implies that on the domains E\ respectively Ef (recall that Ef repre- 
sents Ei in the adapted coordinates y) the conditions y\ ~ 2~\ y 2 ~ 2~ k are equivalent 
to the conditions 

x x ~ 2~ j and <f>(x) ~ 2~ {A ' j+Blk) 

(cf. Lemma 6.1 in [41]). 

Working in the coordinates y, after re-scaling of the measures fij^ to get normalized 
measures Vjj. supported on a surface Sj 7 k where y\ ~ 1 ~ y 2 , by means of the formula 
above one eventually finds that Sj t k is a small perturbation of the limiting surface 

Soo := {(yi, yT^iO), cyf'y 2 l ) : yi ~ 1 ~ 2/2}, 

But \d{cyi l y 2 l ) / dy 2 \ ~ 1, since B\ > 1, which shows that S^, and hence also Sj^, is a 
smooth hypersurface with one non- vanishing principal curvature (with respect to y\) 
of size ~ 1. This implies that 

<c(i+iei)- 1/2 , 

uniformly in j and k. Applying Greenleaf 's restriction theorem to these measures, and 
scaling these estimates back, we eventually arrive (in a not completely trivial way) at 
the estimates (5.8). It is important to observe here that Greenleaf s results implies 
restriction theorems for p' > 2(1 + 2) = 6, which is sufficient for our purposes, since 
p' c >2d + 2, where d > 2. 
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Contribution by the domains D\. Let us next turn to the domains D[. After dyadic 
decomposition of the domain Di in the adapted coordinates y by means of the k 1 - 
dilations and suitable re-scaling, the re-scaled measure V\. corresponding to the mea- 
sures fik turns out to be of the form 

and by means of a finite partition of unity, we may assume that the amplitude fj is 
supported in a sufficiently thin set U(c ), on which 

yi ~ 1 and \y 2 - c Q y{ l \ < ey\ l . 

This measure v k is supported in a variety Sk which in the limit as k — > oo tends to the 
variety 

Soo := {<7oo (2/1,2/2) := w(0)v? ■ (2/1,2/2) e Ufa)}, 

since mK,[ — K 2 < aiK\ — K l 2 = and since fc tends to Here, Co is fixed with |co| < iVj. 

Again, we have to prove uniform restriction estimates for the family of measures z/^. 
Depending on c , different cases may arise. 

1. Case. 0*2^,(1, Co) 7^ 0. Then we may use z 2 := <%(2/i, 2/2) in place of y 2 as a new 
coordinate for (which thus is a hypersurface) , and since y\ ~ 1 on C/(cq), we find 
that Sqo, hence also S^, is a hypersurface with one non-vanishing principal curvature. 
Then we may essentially argue as for the domains E\. 

2. Case. <9 2 0",(l,c o ) = 0, but <9i0^(l,co) 7^ 0. In this case, since is a k z - 
homogenous polynomial, by Euler's homogeneity relation we have also Co) 7^ 0. 
One can then show that one can fibre the variety S^ into the family of curves 

7e(2/i) := £00(2/1, q/?) = (2/i,^(O)j/™,0" i (?/i,c?/i i )), 

1/k' 

for c sufficiently close to Co, and one finds that the curve 70(2/1) = (2/1, w (0)2/i\ &o2/i *) 
has non-vanishing torsion, since 60 7^ 0. The same applies then to the curves j c , and 
for k sufficiently large, we do obtain the analogous results for the varieties Sk- 

This allows to decompose the measure dv^ as a direct integral of measures dT c 
supported on curves % with non-vanishing torsion. 

We may thus apply Drury's Fourier restriction theorem for curves with non- vanishing 
torsion (cf. Theorem 2 in [23] and [5], [20]) to the measures dT c and obtain uniform 
estimates 

(| |/| 2 dT c )*< CpII/IIlpcrs), 

when p' > 7 and 2 < p'/Q. Since we assume here that p' c > 2(d+l) > 2(5+1) = 12, these 
estimates, after re-scaling to the measures /^~, yield the desired restriction estimates 
for the contributions by the domains D\. 

Notice that it is here that we need the condition d = hn n > 5. 
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3. Case. d20 a i(l,co) = and <9i0" ; (l,co) = 0. Then Euler's homogeneity relation 
implies that also <%(1, Co) = 0, so that (p a K i has a real root of multiplicity B > 2 at 
(1, c ), and one finds that 

(5.9) rAvu V2) = vf{v2 - c y?) B Q(yi,y2), 

where Q is a ^'-homogenous smooth function such that Q(1,cq) ^ and Q(1,0) ^ 0. 
One can also prove that B < d/2. 

We can then essentially follow the Stein- Tomas method for proving L p -L 2 - restriction 
estimates. We localize to frequencies of size A > 1 by putting 



where xi is a smooth bump function supported where |£| ~ 1. We claim that the mea- 
sures v£ satisfy the following estimates, uniformly in k > k , provided k is sufficiently 
large and e' sufficiently small: 

(5.10) Halloo < CA- 1/B ; 

(5.H) KIU < CA 2 ~ 1/B . 

Indeed, 



^2< m 4-4)'=j /2 +^a;(2- K iSi)) +&<t>k(v) 



.A 

which, in the limit as k — > oo, simplifies as 

2(0 = Xi (|) / e^"«rtW«,W)] m dy . 

Now, if |6| > c|(£i,£2)|; then an application of van der Corput's lemma to the inte- 
gration in y 2 yields |z^(OI ^ |6|~ 1/B (cf. (5.9)), and if |6[< 1(6,6)1, we may apply 
van der Corput's lemma to the ^-integration and obtain < |(6,6)l -1 ^ 2 . Since 

B > 2, and because van der Corput's estimates are stable under small perturbations, 
we thus obtain (5.10). 

In order to verify (5.11), observe that u^(xi, x 2 , X3) is given by 



A J (J 7 Xi)(A(x! - yi),A(x 2 - u(0)y?), A(x 3 - <t>° t ( yi , y 2 ))fj(y) d Vl dy 2 , 

hence, by a change of coordinates, 

1^(^1,^2,^3)1 < A 2 J p(zi) p(A(x 3 - (f) a K i{xi - ^,2/2)) 771(^1 - ^y 2 )dz x dy 2 , 

where p and r)\ are suitable, non-negative Schwartz functions, and r\\ localizes again to 
U(co). However, since \d 2 4>^i(yi,y 2 ))\ ~ 1 on the domain of integration, classical sub- 
level estimates, originating in work by van der Corput [16] (see also [1], and [13], [26]), 
essentially imply that the integral with respect to y 2 can be estimated by 0(A -1 / B ), 
uniformly in y\ and A. 
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Interpolating the estimates (5.10) and (5.11), and applying the standard Stein- Tomas 
argument (see, for instance, [34]), it is easily seen that we can sum the corresponding 
estimates over all dyadic A ^> 1, and we obtain the L p -L 2 restriction estimate 

( J \f?dv k ) 1/2 <C p \\f\\ LP 

whenever p' > AB, uniformly in k, for k sufficiently large. Since B < d/2, we have 
p' c > 2d + 2 > 4B, so that the range p > 4B does include the critical value p = p c . 

Scaling back to the measures /i^, we find that also the original measures fik satisfy a 
uniform restriction estimate 

( J \f\ 2 dfi k ) 1/2 <C p \\f\\ LP , 

where C p does not depend on k, provided p' > 2h t + 2. However, this applies to p c , 
since h r ((f>) > hi. 

Finally, observe that we can achieve our dyadic decomposition into the measures \x k 
by means of a dyadic decomposition in the variable xi, so that these uniform estimates 
allow to sum over all k by means of Littlewood-Paley theory applied to variable x\\ 

What remains to be understood is the contribution by the domain _D pr = D\ given 

by 

\x 2 - ^{xi)\ < N x xl x . 

In this domain, the upper bound B < d/2 for the multiplicity B of real roots will 
in general no longer be valid, as examples show, not even the weaker condition B < 
h r ((p)/2, which would still suffice for the previous argument. 

In order to resolve this problem, we apply again a further domain decomposition by 
means of a stopping time argument. 

Notice first that if we are to proceed as for the case I < A, a major problem arises 
only in Case 3 where V0p r (l,co) = 0; in all other cases we can essentially argue as 
before. 

We therefore devise the stopping time argument essentially as follows: 

We put 4>^ := (J) a . If Case 3 does not appear for any choice of Co, then we stop our 
algorithm with (jr- 1 ', and are done. Otherwise, if Case 3 applies to Co, so that coy1 x is 
a root of 0" A , then we define new coordinates z in place of y by putting 

z\ := x\ and Zi := Xi — ip(xi) — Cox" A , 

and express by (p^ in the coordinates z. Again, if Case 3 does not appear (for cj)^ 
in place of (jr- 1 ') in the corresponding z-domain, we stop our algorithm. We also stop 
the algorithm when no further fine splitting of roots does occur, i.e., when there is a 
root with leading term ip{x\) +cqx1 x of which has the same multiplicity as the trivial 
root z 2 = of the principal part of (p^. Otherwise, we continue in an analogous way. 

This algorithm will stop after a finite number of steps, and eventually leads to a 
further domain decomposition of D pr into "homogeneous" domains Dm and transition 
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domains E^, which can eventually be treated by methods similar to those applied for 
the domains Ei and D[. 

5.3. Sufficiency of the condition p' > 2h r (4>) + 2: II. Brief sketch of some ideas 
of the proof when hn n (<j)) < 2. This case turns out to be by far more difficult than 
the case where hn n (<f)) > 5, and its discussion requires numerous further methods and 
ideas which I can sketch only very briefly. 

A first observation is that h r ((j)) = d when d = hn n ((p) < 2, so that 



The starting point of our analysis is the following local normal form of our function 
0. It is closely related to the classification of singularities in [2] and [24] . 

Theorem 5.8. If hn Q ((j)) < 2, then locally is of the form 



Here 6, 6 and ip are smooth, and ip is again the principal root jet, and either 

(a) 6(0,0) 7^ 0, and either bo is flat (singularity of type Aoo), or of finite type n, 
i.e., 6o(xi) = Xi/3(xi), where (3(0) 7^ (singularity of type A n _i); 

or 

(b) 6(0, 0) = and b(xx, x 2 ) = Xibi(xi, x 2 ) + xlb 2 (x 2 ), with 61 (0, 0) 7^ 
(singularity of type D). 

Let us just consider the case where is of finite type A n _i (type D can be treated 
in a rather similar way). 

In a first step, by making use of these normal forms in order to estimate certain 
two-dimensional oscillatory integrals that arise in estimating the Fourier transforms of 
surface carried measures, we can again reduce to the domain (3.7), where \x 2 — cx™| < 
ex™. 

In a second step, if k denotes again the principal weight associated to the principal 
edge of Af(4>), we again perform a dyadic decomposition of the measure \i over the 
domain (3.7), and re-scale in a suitable way. This leads to a phase function 



where 5 = (5 Q , 5\, 5 2 ) = (2^( nKl_1 ) fe , 2 _Klfc , 2~ K2fc ) are small parameters which tend to 
as k tends to infinity, and where b(5iXi, 82X2) ~ 6(0, 0) 7^ and 0(0) 7^ 0. 
What we then need to prove is the following 

Proposition 5.9. Given any point v = (vi,v 2 ) such thatv\ ~ 1 andv2 = v™oj(0), there 
exists a neighborhood V of v in (M+) 2 such that for every cut-off function n G T>(V), 
the measure v$ given by 



p c = 2d + 2. 



(5.12) 



(f>(x 1 ,x 2 ) = b(xx,x 2 )(x 2 - i)(xi)) 2 + 6 (xi). 



(j>(x, 5) := b(5xxi, 5 2 x 2 ) ( x 2 - x™u(6ixi) ) + 5 x±(3(6iXi), 
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satisfies a restriction estimate 

n \T\ 2 du s y /2 <c p jf\\ LP 

whenever p' > 2d + 2, provided 5 is sufficiently small. 

In oder to prove this proposition, we again perform a dyadic decomposition, this 
time with respect to the a^-variable. By means of Littlewood-Paley theory, it then 
turns out that it is sufficient to prove uniform restriction estimates for the following 
family of measures 



{vs,j, f) ■= j f(x, <fi(x, 5)) x(2 2j 0(x, 5))rj(x 1 , x 2 ) dx, 

of the form 

— \ V 2 

\f\ 2 dv 6 A < C p , v ll/H lp(R3)- 

Here, x( x s) is supported on a set where |x3| ~ 1. If 2 2 - ? <5 1, then it turns out that 
this localization means in fact again a localization to a curved rectangle where 

\x\ — V\\ < e and \x 2 — x™u(5iXi)\ ~ 2~ J , 

but in other cases, it has another meaning. 

It turns out that these estimates require a refined spectral decomposition of the mea- 
sures ugjj namely a dyadic decomposition in every dual coordinate ^1,^2,^3- Slightly 
cheating, for every triple A = (Ax,A2,As) of dyadic numbers Aj = 2~ fc % we therefore 
define the functions 1A by 



and choose Xi( s ) supported where \s\ ~ 1 so that ugj = Y^K v fi "where summation is 
essentially over all these dyadic triples A. We have here suppressed the dependency of 

Uj on the small parameters 5. Note that |^| ~ Aj on the support of v^. 

For every fixed A, we then essentially follow again the Stein- Tomas approach, by 

estimating H^Hoo and ||z^||oo- 

This requires the distinction of various cases, depending on the relative sizes of Ai, A2 
and A3. In the end, it turns out that the most difficult case is where Ai ~ A2 ~ A3 and 
2 2j 5o ~ 1, and the main result to be proven in this case is the following. 

Proposition 5.10. Let <fi be of type A n _\, with m = 2 and finite n > 5. Then 



2 

L 14/ll 



(5.13) J2 \f?dvf<C2$ 

2<Ai-A 2 ~A 3 <2 6 J S 

for all j G N sufficiently big, say j > j , where the constant C does neither depend on 
5, nor on j. 
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The proof of this result requires yet further refinements. Indeed, the Fourier trans- 
form of i/j is an oscillatory integral with complete phase of the form 

where a := 2 2 - ? <5 ~ 1 and \$(y,6,j)\ ~ 1. Notice that if |£i| ~ |£ 2 | ~ l&l, then may 
have degenerate critical points, with non-vanishing third derivatives, with respect to 
the variable y\ , so that we encounter oscillatory integrals of "Airy type" . 

This case requires a further dyadic frequency decomposition with respect to the dis- 
tance to certain "Airy cones." For the quite comprehensive details, I refer to [41]. 

The case where 2 < d < 5 turns out to be the most delicate one, since on the one 
hand, we can no longer apply Drury's restriction theorem, which necessitates an even 
more refined domain decomposition, and furthermore, somewhat subtle interpolation 
arguments are needed in many situations in order to handle the endpoint p = p c (cf. 
[42]). 
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